An extended ode solver for sensitivity calculations

A rcccntly published semi-implicit Rung+Kulta algorithm has been extended to compute parameter sensitivities in dynamical systems. The decomposed direct method (DDM) of sensitivity analysis results in a general and compact FORTRAN program, well-suited even for small computers. The straightforward DDM provides more accurate sensitivities than does the numerically efficient, but more involved Green’s function (GF) approach. Results show that for kinetics problems 01” moderate complexity the proposed procedure is competitive with the well known GF method in terms of computer-time requirements.